DETERMINATION OF THE PITCH-ANGLE DISTRIBUTION AND TRANSVERSE 
ANISOTROPY OF INTERPLANETARY PARTICLES 


C. K. Ng 

Dept of Mathematics, University of Malaya 
Kuala Lumpur, Malaysia 


1. INTRODUCTION We present a method to determine the directional 
differential intensity (d.d.i.), expressed in terms of spherical 
harmonics, from sectored particle data, concurrent interplanetary 
magnetic field (IMF) and solar-wind velocity. In Section 2, we show the 
relation between the d.d.i. and the mean sector count rates Xi. In 
Section 3, we show how to estimate the d.d.i. from the measured X* and 
the associated errors due to Poisson statistics. In Section 4, using the 
above method, we determine the pitch-angle distribution and" the 
transverse anisotropy of the d.d.i of low energy protons for the Day 118, 
1978 solar particle event. In Section 5, we discuss an interesting 
correlation between the transverse anisotropy and the IMF direction! 

2. RELATION BETWEEN DIRECTIONAL INTENSITY AND SECTOR RATES 

We express the particle directional differential intensity as 

OO n 

j(p,9,0> = Z Z PK<cos 0) { A„ m (p'> cos rn0 + Br, m (p) sin m0 1 , (1) 

where PS = the associated Legendre functions, p = particle momentum, 
9 - pitch-angle and 0 = gyrophase in the standard coordinate system (Ng, 
1985). When the telescope points in the direction Of,rj), it measures the 
differential count rates 

U(p,Y,T]) = j^j(p,9,0) S [ 6’ (T,rj,6,0) ] sin 0 d0 d0, (2) 

where, following Sentman and Baker (unpublished manuscript), we express 
the angular response function of the particle telescope as 

00 OO 

S(O') = X S k P k (cos ©') = z S k { P k (cos r )P k (COS 0) + 

k “°_ k k -° 

+ 2 ZC(k-m)'/(k+m)*] Pj?(cos K )Pjj(cos 9) cos mCq-0)} . (3) 

It follows from the orthogonality of the spherical harmonics that 

OO „ 

C(p,*,r|) = 4 tvZ [s„/(2n+l)] Z{ Ar»m (p)cos mq + B„ m (p)sin mqj Pf?(cosT ).(4) 

n-O m-O^ 

As the telescope sweeps over sector i, we average (4) to obtain X 4 (p) the 
mean differential count rate over sector i, 

°° r " 

X»(p) = 4 7T Z [S„/(2n+l)] z { A„ m (p)< Pr(cos Y)cos mri > t + 

n-O m-O 1 

+ B nm (pX PR(cosV)sin mq > t } , (5) 

where < >i = ( 1 / Ayr) + 1 dp, 'p = the spacecraft (s/c) spin-angle 

measured from the y <- projection of the IMF onto the spin plane, and 
( V^i >y r i+i) defines sector l (see the 2nd coordinate system in Fig. 1). 
For multiple-telescope systems (Sanderson & Hynds, 1977), eqn (5) should 
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be repeated for each telescope. 

We now illustrate by specialising ( 1 ) and (5) to 8-sectored data for 
a telescope sweeping m the spacecraft's spin plane: 

* a 

j(p,0,0) = X A„©P„(COS 0) + X A„iPA(COS 0) COS 0 f (6) 

4 n-o n_1 3 

X* = 4tt X S r ,A„o<3?,o/(2n+l) + 4 Tr X S„A nl 0‘ x ! (2n+l ) , <7) 

mi— O n— 1 

n 

where Q^o = < X a” sin J 0 B cos J y>>i, (8) 

J-O ^ 

OAi = <(W xl cos Y' ,+ Kl 2 sin yd X b3 sin J 0* cos J yr> t , (9) 

J-o 

a3 = coefficient of x J in P„(x), b3 = coefficient^ x^ in ^PAfxl/d-x 2 ^ 

0 B = angle between IMF and s/c spin-axis, and (W x i, W X2 , W X3 ) denotes a 

unit vector in the direction of ExB in the 2nd coordinate system 
(Fig. 1). Note that the integrations in (8> and (9) may be performed 
readily in closed form. 

3. DETERMINATION OF j FROM MEASURED X* To simplify notation in the 


following, 

let 

D„ = A„ 0 (n=0,4) , 

R‘ = Q, 1 ,© (n=0,4) , 

Dm 55 A>v- 4 # 1 
Rr> “ Om — -4 , 1 

(n=5, 7) , 
(n=5, 7) . 


(10) 

(11) 

7’ 

We least-square fit X± = Z C n 

r\ — O 

Ri, (i=0,7), 



(12) 

to the 8 measured sector rates Xt. 

This yields 





7 7 

Z HnjCj = X RnX* 

J-O i-O 

f 

(n=0, 7) , 


(13) 

where 

H mJ = XRJ.Ri , 
1 

(n=0, 7 

; j=0, 7) . 


(14) 

Inverting 

(13), we have 

c„ = Z 

Z RfcX* , 

i 

(n=0, 7) , 


(15) 

and thence 

D„ = X Mnl 
1 

Xi, 

(n=0, 7) , 


(16) 

where M„i 
Oil , and 

is ultimately expressed in terms of 
Sr,. Using (16), we may determine 

QriOf 

the 

3 

I y 



coefficients Ar, m in (6) by a matrix multiplication 
into the measured sector rates. When the IMF pro- I 

jection lies on a sector boundary, the symmetric 
matrix H„j becomes singular. So we drop the A.*© 
term in (6) whenever the IMF projection comes 
within 2° of a sector boundary. 

Suppose that K* counts are registered over the 1 
time interval t c in sector l. Assuming Poisson 
statistics, we estimate 

X* = K t /tc, var(Xi) = Ki/t?, var(D„) = XMSivar(Xi), 


I s 


Fig. 1 


If we define the anisotropy - D„/D©, then, providing the counts are 
not too low, _ _ 

var<£„) - ( 1/D§) X (M„i - Moi) a var(Xi). ( 18 ) 


Systematic errors are, of course, much harder to estimate. 
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D Ay 119, 1978 

2i 15-mmute intensity and anisotropies of 1.4-2. 5 MeV protons 

APPLICATION As an example, we show in Fig. 2 the 15-min averages 

£>io, ^aor ^ao( and -£»» for 1.4 - 2.5 MeV protons, determined 


• Of 
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with the above method using 8 -sector ed particle data (P.I.: R.E. Vogt, 
CalTech), concurrent IMF (P.I.s N.F. Ness, 6 SFO, and hourly solar wind 
speed (P.I.s H. Bridge, MIT), measured aboard IMP -8 and accessed through 
NSSDC. Some typical standard errors due to Poisson statistics only are 
indicated by vertical bars. 

Estimating the spectral slopes by using the corresponding results 
for 4-13 MeV protons, *we have found the Compton-Gett ing corrections <Ng, 
1985) for transformation to the ExB drift frame to be small, (.002 A c .o, 
<.01, <.01, <.02 for Aoo, *io, and £30 respectively. (For transfor- 
mation to the solar wind frame, the corrections are of the order of 
0.04 Aoo* 0.1, 0.3, and 0.5 respect i vely) . Thus Aoo* %xo f \zo f essen- 
tially characterise the pitch-angle distribution in the ExB drift frame. 

5. DISCUSSION The Compton-Gett ing correction, tWj.(3 - pA&o/Aoo + 
pA£o/5Aoa) to the transverse anisotropy -£n is indicated by the dots in 
Fig. 2. The observed anisotropy varies in phase with this correction but 
is much larger iri magnitude. The same feature, even more marked, is seen 
for 4-13 MeV protons. What is the cause of this large discrepancy*? 

For Fig. 2, IMP- 8 's GSE coordinates in Re varies from (21.6, 21.5, 
5.5) to (5.7, 29.5, 18.1). The times when the IMF is connected to the 
nominal bow shock (BS) are indicated by horizontal bars in Fig. 2 (Ng & 
Roelof, 1977). At ~1537 UT Day 119, some solar particles with guiding 
centres below the IMF through IMP -8 are probably shadowed by the nose of 
the BS, resulting in the observed peak value of -^n = 0.53. However BS 
connection does not account for the general variation of -£>n in Fig. 2. 

Close inspection reveals that the sector plot of Xi lags behind the 
IMF in directional changes. Hence a tentative explanation is that some 
observed 15-rnin averages contain a substantial fraction of non-gyrotropic 
distributions which reside a short distance (~1 gyroradius) beyond a 
'kink' in the IMF. An alternative explanation is as follows. When 0» 
swings rapidly in an averaging interval such that its average value is 
close to one end of the range of values, then a field-aligned anisoptropy 
£io would "induce" a nonzero value of -^n, i.e., the apparent value 
of -£,11 is not real. Further studies with smaller averaging intervals 
would clarify this matter. 

6. CONCLUSION We have shown how to obtain the directional differen- 
tial intensity referred to the standard coordinate system (Ng, 1985) from 
sectored particle data and concurrent IMF and solar wind data. The 
corrections for transformat ion to the ExB drift frame are explicitly 
calculated and found to be small for these -1.5 MeV protons. However, 
the new correction formulae would be important for < 500 KeV protons. It 
is tentatively suggested that the ' observed ' transverse anisotropy may in 
large part be induced by a rapidly changing IMF m the presence of a 
f l el d-al i gned am sot r opy . 
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